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Parameterization of the antiproton inclusive production cross section on nuclei 
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A new parameterization of the p inclusive production cross section in proton-proton and proton- 
nucleus collisions is proposed. A sample of consistent pA — > pX experimental data sets measured on 
1 < A < 208 nuclei, from 12 GeV up to 400 GeV incident energy, have been used to constrain the 
parameters. A broader energy domain is covered for the pp — > pX reaction and with a simplified 
functional form used in the fits. The agreement obtained with the data is good. The results are 
discussed. 
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PACS numbers: 13.85.Ni Inclusive production with identified hadrons 



I. INTRODUCTION 

An accurate description of the inclusive antiproton 
production cross section in proton-nucleus collisions, nec- 
essarily relies on the empirical approach to the experi- 
mental data since theoretical calculations can provide at 
best approximate values in the current stage of the the- 
ory. The aim of the present work was to develop a handy 
analytical parameterization for the description of the in- 
clusive p production cross section in proton-proton {pp) 
and proton-nucleus (pA) collisions on the basis of the ex- 
isting body of data, updating the former works on the 
subject. 

The motivations of the work have their origin in the 
needs of Cosmic Ray (CR) physics where a good knowl- 
edge and a good description of the p inclusive production 
cross section is a key requirement for a detailed under- 
standing of the production and propagation of secondary 
galactic and atmospheric antiprotons. The p component 
of the CR flux is an important window for cosmology. 
The main contribution to this flux originates from the 
interaction of the high energy CR flux with the interstel- 
lar matter (ISM) in the galaxy. These ps are called sec- 
ondary galactic. In addition to the secondary products, a 
primary component could exist, undergone for example, 
by the dark matter neutralino annihilation 0] or by the 
evaporation products of primordial black holes , both 
being of major physical and astrophysical interest. Such 
signatures could be obtained only if the basic processes 
of galactic and atmospheric p production cross section in 
pp and pA collisions are known with a good enough ac- 
curacy over a momentum range extending from around 
threshold up to a few hundreds of GeV where the CR 
flux becomes vanishingly small. 

The approach used here closely follows the forms used 
by Kalinovskii, Mokhov and Nikitin [3| - referred to 
as KMN in the following - for the description of the 



pA — > pX cross section. The functional form used in this 
reference has been modified in order to reproduce a much 
larger sample of data sets over a much larger dynamical 
range and for a larger range of nuclear mass, than in the 
original work. This work extends a previous effort cov- 
ering a more limited_ domain of incident momentum and 
of nuclear mass 



II. 



INCLUSIVE CROSS SECTIONS IN HADRON 
COLLISIONS 



In high energy hadron collisions the final state is often 
complex, many particles being produced in the collision. 
The inclusive single particle production cross section is 
a quantity of interest in many physics studies, for a re- 
action ab —* cX , where c is the particle of interest and 
where X represents all the other particles produced in 
any quantum final state allowed in the collision. The 
invariant inclusive single particle distribution is denned 
by: 
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where d 3 a/dp 3 is the triple differential cross section for 
detecting particle c within the phase-space volume ele- 
ment d 3 p c . E c is the total energy of c, while pn and p± are 
the longitudinal and transverse components of p c , respec- 
tively. The rapidity variable y — 0.bhi((E+p\\)/(E—p\\)) 
is often used to describe the p» dependence of the cross 
section because of its interesting properties in Lorentz 
transformations 6] . To obtain two last expression in Q , 
an azimuthal symmetry of the differentiel cross section 
was used. It is also convenient to introduce the following 
dimensionless variables: 
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where Xf is the Feynman scaling variable and xr the ra- 
dial scaling variable (which depends only on the radial 
distance from the kinematic boundary jjj), with p* and 
Pfmax being the longitudinal momentum of the particle 
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and its maximum possible value in the center of mass 
(cm) frame respectively, while similarly E* and E^ x are 
the total energy of the inclusive particle and its maxi- 
mum possible value in the cm frame respectively. The 
latter can written as E^ ax = (s — M x in + mi) /2-y/s, 
with M x min = 2m p + mjij 'A being the minimum possi- 
ble mass of the recoiling particle in the considered process 
and i/s the total cm energy. Note that for any nuclear 
collision, the kinematical variables used here will always 
be expressed in the nucleon-nucleon (NN) rather than in 
the nucleon-nucleus cm frame, since the NN cm frame 
is the relevant physical system, the incident nucleon en- 
ergies being on the scale of 10 GeV while the average 
binding energy of the nucleon in the nucleus is about 8 
MeV. Bound nucleons can be considered as free particles 
for the incident protons. 

A parameterization of the inclusive production cross 
section can be guided by some general phenomenological 
features of hadron collisions (See [H, |!j for the general 
Physics context). 

• All experimental hadronic production cross sec- 
tions show a strong exponential decrease in trans- 
verse momentum, the exponential slope being more 
or less incident energy and recoil mass M x depen- 
dent. 

• Hadronic scaling: The inclusive distribution 
/ (ab — ► cX) of particle c, is to a good approxima- 
tion, a function only on p± and Xf (or xn) at the 
high energy limit yfs — > oo. Furthermore, a large 
number of slow particles is produced (low x / val- 
ues), the distribution decreasing rapidly to zero as 
Xf — ► 1, like (1 — xt) . This form can be explained 
by the counting rules in parton model. 

These features are predicted qualitatively by the 
Regge poles phenomenology and the parton model. 

• For a given p±, the inclusive distribution of pro- 
duced particles is (to a good approximation for 
pA collisions) symmetric in the rapidity space with 
f(pj_,y-y cm ) = f{p±,-y + ycm), where y cm is the 
rapidity of the cm in the laboratory frame (Lab), 
or in the cm frame f(p±,y*) = f(p±,—y*). A 
forward-backward symmetry of the cross section 
is expected from first principles (symmetry of the 
wave function and of the interaction) for the NN 
system (In pA collisions however, the absorption of 
low energy particles in the nuclear medium may dis- 
tort the natural symmetry). In the central region 
where y* = 0, the inclusive distribution consists 
of a plateau which width increases slowly with the 
incident energy. This plateau reduces to a simple 
maximum over the energy range considered here. 
The inclusive distribution rises again in the frag- 
mentation region where y* — > ±2/* nax for particles 
which can be produced diffractively but it is not 
the case for p and the inclusive distribution falls 
simply in the fragmentation region. 



III. PARAMETERIZATION OF THE 
p + A -> p + X CROSS SECTION 

Following the approach proposed in Q, the former 
phenomenological features of hadron collisions have been 
used to constrain the parameters of a functional form de- 
scribing the inclusive p production cross section, which 
could reproduce all the relevant experimental data avail- 
able from pp and pA collisions. The data used are listed 
in Table [I] The measurements on nuclear targets cover 
basically the whole range of nuclear mass, from proton 
to lead, over a range of incident energies from 12 GeV up 
to 400 GeV, matching the useful range for CR studies. 

The KMN parameterization used previously Q is in 
very poor agreement with the data listed in Table [H and 
a reexamination of the analytical approach, better con- 
strained by recent data was necessary. The larger in- 
cident energy domain used here required some energy 
dependence to be introduced in the parameterization fol- 
lowing the general features described above as (loose) 
guidelines. 

In this study, the p inclusive cross section will be ex- 
pressed as a function of the three variables, y/s, p± and 
xr (see for example for the relevance of the choice of 
these variables): 
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The following functional form used to describe the p pro- 
duction cross section is an evolved version of the KMN 
formula: 

E ^ = a in A Cil <£)^(l-x R fM^) e -C4*n 
ap d 

C 5 (^f B e- c ^ +C 8 (Vi) C9 e- Clop i] (4) 

where A is the target mass. The total inelastic cross 
section <7, n for pA collisions was borrowed from [T(il | : 



a m {mb) = a [l-0.62exp(-£ mc /200) 

sin(l0.9iC c 28 )] 
a {mb) = A5A - 7 [1 + 0.016 sin (5.3 — 2.63 In A)] (5) 



where Ei nc is the incident kinetic energy in MeV. 

The 10 parameters C\-C\o have been fitted to the set 
of experimental data listed in Table [I] by a standard \ 2 
minimization procedure using the code MINUIT (TT| . 

In relation 10} , the term (1 — xr) C3 originates from the 
hadronic scaling properties, namely, the quark counting 
rules of the parton model of hadronic interactions |8j] (see 
Sect. [n}. It was found empirically in this study, that a 
significantly better result is obtain if the exponent is en- 
ergy dependent. The ln(-y/s) factor multiplying the C3 
coefficient was found to give the best result. The term 
e -C&x R j g m( juced by the Regge regime 0. The last fac- 
tor of relation Q accounts for the transverse momentum 
dependence of the cross section (see Sect.[n}. The anal- 
ysis of the experimental data (Table [1} showed that the 
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term of angular dependence e~ loP± is dominant at low 
energy, E l £ h as 10 GeV, while the term e~ C7P± dominates 
at high energies, E l p ab > 100 GeV. The y/s dependence 
has been introduced to allow the transition from the p± 
to the p\ dependence. The target mass dependence was 

accounted for by the factor A Cl \^) pJ ~, with an en- 
ergy dependent exponent introduced for the same reason 
as above, at variance with the constant exponent used in 
the KMN parameterization. The energy dependence used 
(linear in y/s) accounts for the experimental increase of 
this coefficient with the incident energy found using the 
KMN approach. For incident energies E l ^ h < 55 GeV, 
this coefficient becomes negative. 




FIG. 1: Experimental p production cross section for p + 
Pb collisions at 23.1 GeV jc and at degree scattering angle 
fl^l compared with the parameterization (JIJ, plotted against 
particle momentum p; a j, (left), and rapidity y (right). In this 
latter case the distribution is symmetric around the value of 
the cm rapidity y cm — 1.9, corresponding to po ~ 3.3 GeV/c 
in the Lab. The fit to the cross section data in the upper p > 
po region of the momentum region also determines the values 
of the cross section for p < po where no data is available. 



to piab ~ 0.65 GeV/c, the validity of the fit extending 
likely significantly below this value. This is an important 
point since experimental p cross section data are usually 
scarce below about 1 GeV, and since the production cross 
section for these low energy particles is very important 
because ps originating from neutralino (dark-matter) an- 
nihilation or from the evaporation of primordial black 
holes 2], arc expected within this energy range. 



ANL-Marmer-Cu-12.3 GeV/c 
X z /point = 2.4795 



E ANL-Marmer-Be-12.3 


GeV/c 


X 2 /point = 


.8833 


T ^^^^ 




- i 








-/ /' i 

" 1 1 / 


, , I , 




1.2 1.4 

p„ (GeV/c) 



JEK-Suqaya-Pb-Cu-AI-C-12 GeV/c l J 
X ! /point = 2. 58571 > 




0.6 D.8 1 1.2 1.4 

p„ (GeV/c) 



BNL-Abbott-AI-14.6 GeV/c 
X 2 /point=1 .0331 

1 T 




1.5 2 2.5 3 
PlIi (GeV/c) 



0.6 0.8 
P, (GeV/c) 



FIG. 2: Experimental data points from |l2l Ihjl. Il4j| compared 
to the best fit calculations. 



In Sect. 2, it was mentioned that one of the features of 
the inclusive distribution is its symmetry in the rapidity 
space. By construction our parameterization (relation 
(@J) satisfies this symmetry property since it depends 
only on y/s, p± and xr. This is illustrated in FigQ which 
shows the fit to the 23.1 GeV/c cross section data versus 
particle rapidity (right) and particle momentum (left). 
The maximum of the cross section corresponds to parti- 
cle production with velocity zero in the cm frame, i.e., 
traveling with the cm velocity in the laboratory. The 
corresponding value of this momentum in the Lab is: 
Po ~ y/m p E p /2. Since the cross section distribution is 
symmetric in the rapidity space, the upper branch of ra- 
pidity with respect to the y cm , completely determines 
the values of cross section along the lower branch. This 
is true as well in the Lab frame where the values of the 
cross section for p > po determine the values below this 
momentum. In the case of the figure this means that a 
fit to the experimental values of the cross section above 
Po ~ 3.3 GeV/c also determines the values of the cross 
section below with the same level of accuracy. Fitting 
the data from about yi a b ~ 2.1 up to yi a b ~ 3.2, also de- 
termines the cross section down to yi a t ~ 0.65, i.e., down 



IV. RESULTS FOR NUCLEAR TARGETS 

The data used in the fit procedure are summarized 
in Table [IJ The fit sample included measurements from 
12 GeV up to 400 GeV incident proton Lab energy on nu- 
clear targets going from Deuterium up to Pb nuclei, and 
for momentum transfers up to 6.91 GeV/c. For pp colli- 
sions, the incident cm energy y/s extended from about 6 
up to 63 GeV. 

The \ 2 per point obtained with the parameterized re- 
lation Q is 5.32 ( Table ITTTf) for 654 experimental points 
(see list in Table QJ. The values of the parameters ob- 
tained in the fit are given in Tablc^together with error. 
The correlation coefficients between the parameters de- 
termined in the search are given in appendix. The results 
for nuclear targets (A>2) are shown on Figs [21 to where 
the data points are compared with the calculated values 
are given in appendix. The values of the x 2 per point of 
each set are given in the figures. In each case, some basic 
informations (Authors, beam energy, target nuclei, and 
X 2 per point obtained for the considered set) are given 
on the figures. In all the figures, the top distribution cor- 
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Experience 


Target 


Pine or ^J~s 


p kinematical range 


Slab 






GeV/c (GeV) 


GeV/c 


mrad 


Y.Sugaya et al 


C, Cu, Al, Pb 


12 


p lab : 1.0-2.5 


89 


KEK-PS 1998 [12] 










G.J.Marmer et al 


Be, Cu 


12.3 


p lab : 0.820, 1.030 


0, 0.17, 10 


ANL 1969 [13] 










T.Abbott et al 


Al 


14.6 


m t : 0-0.3 




BNL 1993 [14] 






y: 1.0-1.6 




J.V.Allaby et al 


p, Be, Al, Cu, Pb 


19.20 


pi ab : 4.5-18.3 


12.5-70 


CERN 1970 [15] 










D.Dekker a/ 


p, Be, Pb 


18.8, 23.1 


piab: 4-12 


0, 100 


CERN 1965 [16] 










T.Eichten ef aZ 


Be, Al, Cu, Pb 


24 


Piab: 4-18 


17-127 


CERN 1972 [17] 










P.I.P.Kalmus e£ al 


Be 


24-26 


piab'- 0.6—2.5 


310 


CERN 1971 [181 










J. M. Snow et al 


Pt 


28.4 


pi ab : 0.606-0.730 


0-0.17 


BNL 1985 [19] 










D.S.Barton et al 


p, C, Cu, Al, Ag, Pb 


100 


pi ab : 30-88 




FNAL 1983 [20j 






px: 0.3, 0.5 




J.R.Johnson et ai 


P 


100, 200, 300 


p ± : 0.25-1.5 




FNAL 1978 [21] 






0.05 < x fl < 1.0 




W.F.Baker aZ 


Be 


200, 300 


pi ab : 23-197 


3.6 


FNAL 1974 [221 










J.W.Cronin et al 


Be, Ti, W 


300 


px: 0.76-6.91 


77 


FNAL 1975 [23] 










DAntreasyan et al 


p, d, Be, Ti, W 


200, 300, 400 


px: 0.77-6.91 


77 


T7M AT 1 HTfl TO /1 1 

£JNAL 1979 |24| 










K.Guettler et al 


P 


23 < yfs < 63 


Px: 0.1-0.3 




CERN 1976 [25] 










P.Capiluppi et al 


P 


23.3 < y^s < 63.7 


Piab: 1.5-10 


80-350 


CERN 1974 [26J 











TABLE I: List of the experimental p production cross section data included in the x 2 minimization procedure, in increasing 
energy order. 



parameter 


Ci 


c 2 


c 3 


c 4 


Cs 


value(error) 


0.16990(4) 


10.28(13) 


2.269(7) 


3.707(27) 


0.009205(2) 


parameter 


c 6 


CV 


Cs 


c 9 


Cio 


value(crror) 


0.4812(14) 


3.3600(2) 


0.063940(73) 


-0.1824(15) 


2.4850(6) 



TABLE II: Values of the parameters of relation obtained from fitting the experimental p production cross sections list in 
Table and the corresponding error following the PDG standard conventions. 



system 


parameterization 


X 2 per point 


pp,pA 


KMN [3] 


80.0 


pp,pA 


this work Q 


5.3 


pp 


Tan and Ng [31] 


28.1 


PP 


this work 16J 


3.6 



TABLE III: Comparison between this work and the other pa- 
rameterizations . 



responds to the measured cross section, while each next 
distribution below has been multiplied by a 10 _1 factor 
with respect to the previous one, for the legibility of the 
figure. 



As it can be seen on the figures, the quality of the fits 
varies from fair to excellent. A poor fit is obtained how- 
ever for the 24-26 GeV data from [l^, the calculations 
underestimating the data by a factor of about 2. Never- 
theless, this set has been kept in the fit procedure since 
its contribution is small and not hardly affects the results 
(which is not the case for the data listed in Table llVjl . To 
the opposite, outstandingly good fits have been obtained 
consistently and simultaneously for the CERN data from 
in the 20-25 GeV incident energy range, and for 
the high energy and large momentum transfer data from 

mm. 

Table ITTT1 compares the values of the x 2 P er point ob- 
tained in the present study, with that obtained using the 
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FIG. 3: Experimental data points from |15| compared to the 
best fit calculations. 
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FIG. 5: Experimental data points from [TH ] compared to the 
best fit calculations. 
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FIG. 4: Experimental data points from [Tfj ITsL ITfll] compared 
to the best fit calculations. 



KMN relation 3] for the same data. The latter is seen 
to be more than one order of magnitude larger than the 
value obtained using Q. This gives the scale of the im- 
provements achieved by the present study on the issue. 

These results demonstrate the ability of the proposed 
parameterization to describe the inclusive p production 
cross section on nuclei over the quoted ranges of inci- 
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FIG. 6: Experimental data points from |2(J IgUgH^] com- 
pared to the best fit calculations. 



dent energy, momentum transfer, and target mass, with 
a good accuracy. 
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Experience 


target 


pine/ Einc 


p kincmatical range 


9lab 






(GeV/c) 


(GeV/c) 


(mrad) 


Yu.M.Antipov et al 


Al 


70 


Plab : 10-60 





IHEP 1971[27] 










V.V.A.Abramov et al 


C, Al, Cu, Sn, Pb 


70 


p ± : 0.99-4.65 


160 


TTim i no a Tool 

IMtF 1984 |28J 










W.Bozzoli et al 


Be, Al, Pb 


200 


p lab : 20-37 





CERN 1978 [29] 










I.G.Bearden et al 


Be, S, Pb 


450 


pi ab : 4-8.5 


37, 131 


CERN 1998 [30] 






p ± : 0.11-1.28 





TABLE IV: Antiproton production cross section data not taken into account in the \ 2 minimization procedure, classified by 
increasing energy. See text for explanations. 



10 

i ;f 

io~' 

icT 

1 cf 
- 
10 

1 0~ 

10 



HEP-Antipov-AI-70 GeV/c 
xVpo1nt = 648.577 
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FIG. 8: Experimental data points from [30|l compared to the 
best fit calculations. Data points not included in the final 
search procedure. 



Data discarded from the selection 



The data listed in Table lIVl were not included in the fit 
sample, because of their obvious inconsistency with the 
other data. This is illustrated on Figs. [7| and [H] where 
they are compared to the best fit calculations obtained 



in the previous step on the selected sample. As it can be 
seen, the difference between data and calculated values 
amounts up to about one order of magnitude. The ratio 
goes from 2 to 10 for the Serpukhov experiments [2?], 
|28| . For [2]]], it is about 5, and consistent with a simple 
normalization problem. 

A larger and more surprising disagreement is found 
with some recent CERN data from NA44 |2jj, in par- 
ticular for the measurements in the small rapidity bin. 
Note also that the parameterization © describes quite 
well the data from [24( obtained on the same targets as 
[3fj| over a wider kinematical region (see Fig.[SJ). 

The p cross section data from [2£j appearing in the 
table, were given in the original works in units of the 
corresponding production cross section measured at 
the same momentum. Although the absolute value could 
be obtained using the known 7r~ cross section, the results 
were considered too inaccurate however, and discarded 
from the selected sample. 



B. Analysis of the pp — > pX data 

This reaction is the dominant contribution to the sec- 
ondary p production induced by Cosmic Rays, since the 
interstellar gas is mainly constituted of hydrogen gas. It 
is thus important to obtain as accurate a description as 
possible for the cross section. 

Considering separately the p + p collision data in Ta- 
ble U the parameterization gives for the best fit a 
value of the \ 2 P er point of 7.08. For the same data, the 
well known parameterization of Tan and Ng |3l| gives 
a value of 28.1. In addition, this latter parameteriza- 
tion is valid only for p± = — 0.8 GeV/c and is not 
able to reproduce the large p± data such as those from 
[13 and [I! where p± = 0.76 - 6.91 GeV/c. Note also 
that Tan and Ng's parameterization contains 8 param- 
eters for > 10 GeV {pi ab > 50 GeV/c) and 17 for 
y/s < 10 GeV. 

However, in the course of the study, it appeared that 
some of the parameters of relation (0J had no incidence 
on the resulting fits. The parameterization Q has then 
been revisited and simplified from some of its parame- 
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parameter 


Di 


D 2 


Da 


Di 


D 5 


D 6 


D 7 


value(error) 


3.4610(20) 


4.340(20) 


0.007855(3) 


0.5121(27) 


3.6620(5) 


0.023070(1) 


3.2540(77) 



TABLE V: Values of the parameters of relation © obtained by fitting the experimental p production cross sections list in 
Table 3 for proton-proton collisions and the corresponding error following the PDG standard conventions. 
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FIG. 9: Experimental data points from pp collisions from 
[l5lH^ . l25ll26| compared to the best fit calculations using the 
two parameterizations. See text. 
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FIG. 10: Same as Fig. H for the data from ^[H, The \ 2 
per point for the first three graph is indicated in graph 3. 



ters irrelevant for this particular reaction and from other 
parameters which turned out to be ineffective in the min- 
imization procedure, resulting in the following functional 
form for pp — » pX inclusive production cross section : 



E- 



dp 3 



'in (1 - x R ) 1 e 
D s (VS) D *e- D * 



-D2XR 



D 6 e 



(6) 



In comparison with relation Q , the dependence with the 
mass of the target has been removed since the only proton 
target is considered in this case. In addition, the energy 
dependent factors y^s in front of D\ and D7 in JIJ have 
been also removed because of their ineffectiveness in the 
minimization procedure. 

The parameters D\ to D7 have been adjusted by the 
same \ 2 minimization procedure as previously [Tl| , to the 
set of experimental data listed in Table |U restricted to pp 
collisions. With formula (JSJ, the \ 2 P er point obtained 
for the best fit is 3.59, for 228 experimental point, instead 
of 7.08 with relation (0J. The values of the fit parameters 
obtained with © are given in Table [V] The correlation 
coefficients between the parameters determined in the 
search are given in appendix. Note that the values of the 
coefficients C3 and D2, C5 and D3, respectively, are of 
the same order of magnitude. This was expected since 
they describe the same physics in the relations JIJ and 
©. 

Figs . l9l and 1101 show the pp — > pX data analyzed, com- 
pared with the best fit results obtained for the whole pp 
and pA data sets from Tabled using relation J3J (dashed 
line) and with those obtained for the pp data only using 
relation I© (solid line). The simplified form © clearly 
provides a significantly better account of the measured 
cross sections, the x 2 value obtained being better by a 
factor of about 2 (about 3.6 against about 7). 

The results obtained in this work have been used in 
the calculations of the d, t and He production from p + p 
and p + A collisions in the atmosphere and in the galaxy 



V. ANTIPROTON MEAN MULTIPLICITY 

In this section, the antiproton mean multiplicity, de- 
fined as 



1 



d 3 p 



(7) 



and depending only on ^fs, has been computed by means 
of relations and © and compared with the experi- 
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mental data in pp collision |35j. Note that the origi- 
nal data from [35| have been corrected from the single- 
diffractive contribution to the total inelastic cross sec- 
tion (Ji n j3|J. The corrected antiproton mean multiplic- 
ity should thus be somewhat smaller than the measured 
values (by ~ 15-20%). 




FIG. 11: Antiproton mean multiplicity distribution in the 
whole phase space, calculated using relation @ (solid line), 
and relation (JSJ (dashed line), compared with experimen- 
tal data [BE!, uncorrected (full circles) and corrected (open 
squares) from diffractive contribution. See text for details. 

The results, shown on Fig. ^] are in good agreement 
with the experimental data. Note that below y/s ~ 15 
GeV/c, the results given by the relation (0J and © be- 
come significantly different (about a factor 2 at the max- 
imum) . As expected, the simplified form JfjJ gives a little 
bit better results since the experimental data are from pp 
collisions. 



VI. CONCLUSION 

The parameterization of the p inclusive production 
cross section on nuclei has been revisited by investigating 
a broad collection of data sets available, covering a large 
dynamical domain of incident energy and of momentum 
transfer, for a broad range of nuclear masses. Good re- 
sults have been obtained but for a small sample of data 



sets inconsistent with the other data. The experimental 
p inclusive production cross sections can be reproduced 
to within a few tens of percent over this range, i.e., for 
incident energies from 12 GeV up to 400 GeV, and for 
target mass 1 < A < 208. Theses results constitute a 
significant improvement with respect to the former KMN 
parameterization, decreasing by a factor of about 15 the 
value of the x 2 P er point obtained using the latter. A 
simplified version of the functional form has been devel- 
oped for pp collisions giving also good results up to very 
high energies, much beyond the range dictated by the 
Cosmic Ray Physics requirements which motivated the 
study. This also constitutes an improvement, consistent 
with the data on nuclei, of the Tan and Ng formula used 
so far as a standard in the calculations. 

The parameterisations (@J and @ are also able to re- 
produce the experimental antiproton mean multiplicity 
measured in pp collision with a good accuracy. 

A point to be emphasized is that because of the sym- 
metry of the cross section in the rapidity space, the fitted 
range in the Lab momentum of the particle, usually mea- 
sured above the cm rapidity, also determines the cross 
section at low momenta, a range of major importance 
for Cosmic Ray antiprotons where accuracy is extremely 
important. 
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VII. APPENDIX 

The symmetrical matrix and 10 give respectively 
the correlation coefficients for the parameters C\— Cm 
and Dx—D'j of equations (0} and 10 respectively. 

In relation the coefficients C\ and C2 appear to 
be strongly correlated (correlation coefficient 0.961), as 
it could be expected from their functional dependence. 
On the contrary, coefficient C5 and Cs are not correlated 
(correlation coefficient 0.232), since effective in different 
energy ranges (see Sect. lIIIjl . The same remarks apply to 
the coefficients Di—Dj. 
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/ 1.000 0.961 
1.000 



0.120 -0.200 



0.131 
1.000 



-0.209 
-0.937 
1.000 



-0.148 
-0.199 
-0.321 
0.282 
1.000 



0.128 
0.157 
0.228 
-0.180 
-0.962 
1.000 



0.086 
-0.098 
-0.049 
-0.042 

0.358 
-0.164 

1.000 



-0.067 
-0.067 
-0.655 

0.834 
-0.110 

0.232 
-0.127 

1.000 



-0.048 
-0.060 
-0.620 
0.784 
-0.128 
0.286 
0.007 
0.979 
1.000 



-0.165 \ 

-0.148 

-0.289 

0.311 

0.239 
-0.202 
-0.028 

0.210 

0.148 

1.000 J 



Correlation coefficients for the parameters D\— Dio given in Tables IH1 



/ 1.000 



-0.933 
1.000 



V 



-0.557 
0.654 
1.000 



604 


0.080 


-0.311 


0.087 \ 


730 


-0.141 


0.435 


-0.135 


979 


0.502 


-0.212 


-0.241 


000 


-0.336 


0.033 


0.228 




1.000 


-0.833 


-0.213 






1.000 


0.282 








1.000 / 



(9) 



Correlation coefficients for the parameters G\—C-r given in Tables IVI 
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